Threshold temperature for pairwise and many-particle thermal entanglement in the 

isotropic Heisenberg model 
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We study the threshold temperature for pairwise thermal entanglement in the spin-1/2 isotropic 
Heisenberg model up to 11 spins and find that the threshold temperature for odd and even number 
of qubits approaches the thermal dynamical limit from below and above, respectively. The threshold 
temperature in the thermodynamical limit is estimated. We investigate the many-particle entangle- 
ment in both ground states and thermal states of the system, and find that the thermal state in the 
four-qubit model is four-particle entangled before a threshold temperature. 
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Entanglement is a nonlocal quantum correlation and 
entangled states constitute indeed a valuable resource in 
quantum information processing |Q. Entanglement in 
systems of interacting spins JPthI as well as in systems of 
indistinguishable particles |7|-|ll[ has been investigated. 
In particular entanglement in both the ground state [^]||] 
and thermal state of a spin-1/2 Heisenberg spin 

chain have been analyzed in the literature. The intrigu- 
ing issue of the relation between entanglement and quan- 
tum phase transitio n |jl2| have been addressed in a few 
quite recent papers [ |13yi4| ]. 

The state we will consider is the thermal equilibrium 
state in a canonical ensemble. The system state at fi- 
nite temperature is given by the Gibb's density operator 
p T = exp {—H/kT) /Z, where Z =tr[exp {—H/kT)] is the 
partition function, H the system Hamiltonian, k is Boltz- 
mann's constant which we henceforth will take equal to 1, 
and T the temperature. As px represents a thermal state, 
the entanglement in the state is called thermal entangle- 
ment Q. It is important to stress that, although the 
central object of statistical physics, the partition func- 
tion, is determined by the eigenvalues of H only, thermal 
entanglement properties generally require in addition the 
knowledge of the energy eigenstates. In some models such 
as the isotropic Heisenberg qubit rings the pairwise ther- 
mal entanglement can be completely determined by the 
partition function fl5|j . 

Let us briefly review the main results of Ref. ]ll| . We 
consider a physical model of a ring of N qubits interacting 
via the isotropic Heisenberg Hamiltonian 



H 



N-l 

J <Sj,i+i + JSn,i, 

i=l 



(1) 



where Sj,j+i =i (1 + c?j ■ is the swap operator be- 

tween qubit i and j, &i = (<7 ia ., <j iy , <j iz ) is the vector of 
Pauli matrices, and J is the exchange constant. In Ref. 
Jj"5| we proved that there is no thermal entanglement be- 
tween two nearest-neighbour qubits in the ferromagnetic 
(J < 0) isotropic Heisenberg rings at any temperature. 
So we only consider the antiferromagnetic case and the 



exchange constant J is then assumed to be 1. A di- 
rect relation is established between the concurrence |l6| ] 
quantifying two-qubit entanglement and a macroscopic 
thermodynamical function, the internal energy U Jig]: 



C(N) = max ( 0, -— ) = max ( 0, — ^ 
1 ' NJ \' NZ d/3 



(2) 



where (3 = 1/T and C(N) refers to the concurrence for 
two nearest-neighbor qubits of the A^-qubit ring. The en- 
tanglement is uniquely determined by the partition func- 
tion of the system. Note that there is a slight difference 
between the expression of the concurrence and the one 
given in Ref. due to a constant term added in Eq. (jl]) 
for convenience of numerical calculations. The reason 
why the entanglement is determined only by the eigen- 
values of the system is that there are many symmetries in 
the isotropic Heisenberg model such as translational in- 
variance and rotational symmetry along x, y, z directions, 
et al. If we apply a nonzero magnetic field along z direc- 
tion we only have rotational symmetry along z direction. 
In this case we can not determine the entanglement by 
eigenvalues alone. 

From Eq.(||) the concurrence for the ground state (GS) 
with energy Eq$ is then followed by 



(3) 



which is obtained in Ref. for even-number qubits. But 
now we know that it is valid for any number of qubits. 
From the available results 17 - 20| of the ground-state en- 
ergy we can determine the concurrence directly. For the 
case of even-number qubits the concurrence is given in 
Ref. H). Here in Table I we give the concurrence of 
ground states for N qubits up to A = 11. From the 
table we find that the entanglement increases with the 
increase of odd number N. On the contrary the entan- 
glement decreases with the increase of even N. We also 
see that the ground state is always entangled except the 
case of N = 3. 
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Since dU/dT > and Cgs(3) = there is no entan- 
glement at any temperature for the case of A = 3. For 
other cases, Ccs(N) > (A ^ 3), and the concurrence 
(|J) is a monotonic decreasing function of the tempera- 
ture. It is then obvious that there exists a A-dependcnt 
threshold temperature T t h(N) after which the thermal 
entanglement disappears. The threshold temperature is 
determined by the equation U(T t h) — 0. Therefore the 
problem left is to examine the threshold temperature and 
we address it in this brief report. We will determine the 
threshold temperature both analytically and numerically. 

The number of qubits we will consider is up to A = 1 1 . 
Analytical results can be obtained for very small rings. 
For A — 2, the eigenvalues of H is easily found to be 
-2(1) and 2 (3) where the number in the parenthesis 
is the degeneracy. For A = 3 the eigenvalues of H are 
0(4) and 3(4) [||. For A = 4, by direct diagonalization, 
the eigenvalues are given by -2(1), (3), 2(7) and 4(5) 
(see later discussions). From the eigenvalues the partition 
functions Z(N) (A = 2,3,4) simply follow as 

Z{2) = e 20 + ier 2 ^ 

Z(3) = 4 + 4e~ 3/3 , 

Z(4) = 3 + e 2/3 + 7e^ + 5e-^. 

The substitution of Eq. (||) into (||) leads to 
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T th (N)\ N odd < T th (oo) < T th (N)\ N 



(6) 



(4) 



C(2) = 
C(3) = 



C(4) = 



Z(2) 
1 

W) 

1 



2/3 



max (0, e 



max (0, — 4e 



3e 



-20\ 



-3(3\ _ 



■ max 



(0, 



7e 



-2/3 



lOe 



-4/3\ 



(5) 



2Z(4) 

which are explicit analytical expressions for the concur- 
rence. The threshold temperatures are directly obtained 



4/ In 3 



3.641, 



from the above equation: T t h{2) = 
T th {3) = 0, and T th (4) = 1.726728. 

For A > 5 we use numerical method to determine the 
threshold temperature. We mention that some analytical 
results of the eigenvalue problem for A = 5, 6 and 7 can 
be found in Ref. Bethe's Ansatz |22| can be used 

for the calculations of eigenvalues. Here we use the ex- 
act diagonalization method to calculate the eigenvalues. 
The obvious symmetry [H, S z ] — allows to decompose 
the Hilbert space of the system into a set of orthogonal 
subspaces. Here S z = is the collective op- 

erator for A qubits. The dimension of the subspace is 
much smaller than the original Hilbert space, and then 
we can diagonalize the Hamiltonian in the subspaces. 

The numerical results for the threshold temperature 
are given in Table II and III, from which we find that 
threshold temperature increases (decreases) as A in- 
creases for the case of odd (even) qubits. We also observe 
that the threshold temperature converges very quickly. 
Then we can estimate the threshold temperature in the 
thermodynamic limit T t h(oo) rs 1. 5903 from the numeri- 
cal results of these finite systems. So we have the follow- 
ing relation 



which implies that the threshold temperatures for A odd 
and A even approach the thermal dynamical limit from 
below and above, respectively. 

We have considered only pairwise thermal entangle- 
ment of a multiqubit system. Another type of entangle- 
ment is the A-particle entanglement which involves all 
A particles. As far as we know there are no discussions 
on many-particle thermal entanglement. In the following 
we will consider the many-particle entanglement in the 
ground states of the Heisenberg models up to 11 qubits 
and then study the many-particle thermal entanglement 
within the four-qubit Heisenberg model. 

For a three-qubit pure state, there exists a good mea- 
sure of three-particle entanglement |^3| . Here we use the 
state preparation fidelity F of a A-qubit state p to inves- 
tigate the many-particle entanglement. The state p can 
be either pure or mixed. The fidelity is defined as pf| 



F(p) = (* G hz|p|*ghz) 



(7) 



where |* G hz) = l/V2(|00...0) + |ll...l)) is the A-particle 
Greenberger-Horne-Zeilinger (GHZ) state p4j . 

The sufficient condition for A-particle entanglement is 
given by p5) 



F(p) > 1/2. 



(8) 



We have a freedom to choose other GHZ states such as 
|*ghz> = 1/V2(|00...01) ± |11. ..10)), etc. By local uni- 
tary operations we can transfer these states to the orig- 
inal GHZ state and the operations do not change the 
entanglement. Detail discussions on this sufficient condi- 
tion can be found in Ref. 
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Let us study the many-particle entanglement in the 
ground states of the Heisenberg models for the number 
of qubits up to 11. It is known that |l9) the ground 
state which has S z — 0, is non-degenerate if N is even, 
and the ground state which has S z = ±1/2, is four- fold 
degenerate if N is odd. Then the ground state is a pure 
(mixed) state for even (odd) N. We numerically find the 
ground state of the system, from which we calculate the 
fidelity. The form of the GHZ state is choose in order 
that the fidelity is as large as possible. The numerical 
results is given in Table IV. First we observe that the 
fidelity of four-qubit ground state is larger than one half, 
hence it is a many-particle entangled state. The fidelity 
is less than one half for N from 2 to 11 except N = 4. If 
we consider even or odd number A, the fidelity decreases 
as A increases. From the analytical expressions of the 
ground states for A = 3 || and A = 4 (see below), we 
can analytically obtain the fidelity F=l/6forA = 3 
and F = 2/3 for A = 4, respectively. 

Having found that the ground state of the four-qubit 
Heisenberg model is a 4-particle entangled state, we ex- 
pect that the thermal state will be also a 4-particle en- 
tangled state for small temperature. We also expect a 
threshold temperature before which the ground state is 
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4-particle entangled. To study the many-particle ther- 
mal entanglement we need to know all the eigenvalues 
and eigenstates. For convenience following notations are 
used 



|n)o = r n - 1 |l) ,|l) = |1000), 

|n>! = 7*- 1 |l) 1 (n = 1,2,3,4), II)! = |1100), 

\n) 2 = T n - 1 \l) 2 (n = 1,2), |1) 2 = |1010). 



(9) 



Here T is a unitary cyclic right shift operator fl9[ | 
defined by its action on the basis T\mim2m3m4) = 
\m4m1m2m3} and \m1m2m3m4) = \mi) £5 \m 2 ) ® Im^) (g> 
|m 4 ). 

A direct diagonalization of the four-qubit Hamiltonian 
gives all eigenvectors and eigenvalues, 

E = 4,|*o) = 10000), 
Em = 4,|*i 5 ) = 11111), 

1 4 

E k = (2 + t k + t,- 1 ), |* fe ) = - £ f-» , 

71=1 

! 4 

s fe+ 4 = (2 + < fc + 1^ 1 ), i* fc+4 ) = - tr^Mo, 



Eg 

E10 
En 

E\2 
E13 

En 



-2,|*9) 

0,|*io) 
2,|*n) 
2,|*i2) 
2,|*i3) 
4,|*i4) 



2-y/3 



^r"- 1 |i)i-2^r- 1 |i) 2 



^=(11010) - 10101)), 
-^(|1100) -10011)), 



7f (|1001) 
5(11100)4 



- 10110)), 

10011) - |1001) - |0110)), 



v \n—l n—1 / 



F(pr) 



e- 4 ^/3 + 2e 2 ^/3 

3 + e 2/3 + 7e -2/3 + 5e -4/3 ' 



(11) 



T t /j Ri 0.83. Since the condition F > 1/2 is a sufficient 
condition we can not say that the thermal state has no 
four-particle entanglement after T t h- From the figure we 
can see that the thermal state has both pairwise entangle- 
ment and four-particle entanglement before the threshold 
temperature. But both entanglement arc not maximal. 

In conclusion we investigated the threshold temper- 
ature for the pairwise thermal entanglement in the 
isotropic Heisenberg model with finite number of qubits 
up to N = 11. We find that the threshold temperatures 
for odd and even number of qubits approach the ther- 
mal dynamical limit from below and above, respectively. 
Although we only make calculations of the threshold tem- 
peratures for small spin rings we can well estimate the 
threshold temperature in the thermodynamical limit as 
Tth(oo) ~ 1. 5903 since T t h{N) converges very quickly as 
N increases. By using one sufficient condition for many- 
particle entanglement, the state preparation fidelity is 
larger than one half, we investigated the many-particle 
entanglement of both ground states and thermal states 
in the Heisenberg models. We find that the ground state 
of the four-qubit model is a four-particle entangled state, 
and also find that there exists a threshold temperature 
Tth ~ 0.83 before which the thermal state is a four- 
particle entangled state. It is interesting to see that pair- 
wise entanglement and four-particle entanglement (not 
' maximal) coexists in this system. 
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(10) 



where the parameters tf. = exp(i7r/c/2)(fc = 1,2,3,4), 
and the operator A^ = erf 4 = a x ® a x ® a x ® a x . 

These eigenvalues and eigenstates completely deter- 
mine the thermal state pr- We let our GHZ state be 
|*GHz) = l/\/2(|1010) + |0101)). Then from Eq.(@) and 
all the eigenvalues and eigenstates we get 



It is easy to verify that F = 2/3 in the limit of T — * as 
we expected. 

In Fig.l we give the numerical result for the state 
preparation fidelity against temperature. For compari- 
son we also plot the concurrence. First we observe that 
there exists a threshold temperature T t h for the state 
preparation fidelity before which the thermal state px is 
four-particle entangled state. The threshold temperature 
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TABLE I. The ground-state concurrence for the number of 
qubits from 2 to 11. 



A 



6 



7 



8 



9 



10 



11 



Cos 1.0 0.0 0.5 0.247 0.434 0.316 0.412 0.344 0.403 0.358 



TABLE II. 
qubits 



The threshold temperature for even-number 



N 



6 



8 



10 



T th 3.64095691 1.72672823 1.60976354 1.59167655 1.59038369 



TABLE IV. The state preparation fidelity of the ground 
states for the number of qubits from 2 to 11. 



N 


Fidelity 


GHZ states 


2 


1.0000 


[01) - 


10) 


3 


0.1667 


010) + 


101) 


4 


0.6667 


|0101) + 


1010) 


5 


0.0873 


|01010) + 


10101} 


6 


0.4580 


010101) - 


101010) 


7 


0.0505 


0101010) + 


1010101) 


8 


0.3173 


01010101) + 


10101010) 


9 


0.0306 


1010101010) + 


101010101) 


10 


0.2205 


0101010101) - 


1010101010) 


11 


0.0191 


101010101010) + 


10101010101) 




FIG. 1. The state preparation of fidelity (solid line) and 
the concurrence (circle dots) against the temperature. 



TABLE III. The threshold temperature for odd-number 
qubits 



N 


3 


5 


7 


9 


11 


Tth 


0.0 


1.53825517 


1.58556286 


1.58979598 


1.59020107 
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